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Abstract 

We prove several theorems concerning Tutte polynomials T{G, x, y) for recursive families 
of graphs. In addition to its interest in mathematics, the Tutte polynomial is equivalent 
to an important function in statistical physics, the Potts model partition function of the 
g-state Potts model, Z{G,q,v), where n is a temperature-dependent variable. We determine 
the structure of the Tutte polynomial for a cyclic clan graph G[(Kr)m, L = jn] comprised 
of a chain of m copies of the complete graph such that the linkage L between each 
successive pair of K^s is a join jn, and r and m are arbitrary. The explicit calculation of 
the case r = 3 (for arbitrary m) is presented. The continuous accumulation set of the zeros 
of Z in the limit m ^ oo is considered. Further, we present calculations of two special 
cases of Tutte polynomials, namely, flow and reliability polynomials, for cyclic clan graphs 
and discuss the respective continuous accumulation sets of their zeros in the limit m —>■ cx). 
Special valuations of Tutte polynomials give enumerations of spanning trees and acyclic 
orientations. Two theorems are presented that determine the number of spanning trees on 
G[{Kr)m, jn] and G[{Kr)m,id], where L = id means that the identity linkage. We report 
calculations of the number of acyclic orientations for strips of the square lattice and use these 
to obtain an improved lower bound on the exponential growth rate of the number of these 
acyclic orientations. 
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1 Introduction 


1.1 Tutte Polynomial 


The Tutte polynomial [Q]-P| (sometimes called the dichromatic or Tutte/Whitney polyno¬ 
mial) contains much information about a graph and includes a number of important functions 
as special cases. Some additional early papers and reviews are |P]-[T^. In this paper we shall 


present a number of results on Tutte polynomials for recursive families of graphs and on re¬ 
lated asymptotic limiting functions in the limit where the number of vertices on these graphs 
goes to inhnity. In this introductory section we give some relevant dehnitions, notation, and 
previous relevant theorems. 


Def. Let G = (V, E) be a graph with vertex and edge sets V and E and let the 
cardinality of these sets be denoted as \V\ and \E\, respectively. One way of dehning the 
Tutte polynomial of a graph G is 


T{G,x,y) = J2 (i- lyo') (1,1,1) 

G'CG 


where G' = {V,E') is a spanning subgraph of G, i.e., E' C E, k{G') denotes the number of 
connected components of G', and c{G') denotes the number of independent circuits (i.e. the 
co-rank) in G', satisfying c(G') = \E'\ +k{G') — |Id|. The hrst term in the summand of (|1.1.1|) 
can equivalently be written (x — ) where the rank r(G) of a graph G is defined as 

r(G) = \V\ — k{G). We shall be interested in connected graphs here, so that k{G) = 1. 

An equivalent way to define the Tutte polynomial is as follows: if El is obtained from a 
tree with \E\ edges by adding i loops, then 

T{H,x,y) = x^%^ ( 1 . 1 . 2 ) 

and if e G is not a loop or a bridge (= isthmus or co-loop), then T{G,x,y) satishes the 
deletion-contraction property 


T(G, X, y) = T{G - e, x, y) + T{G/e, x, y) (1.1.3) 

where G — e = G\e denotes the graph with the edge e deleted and G/e denotes the graph 
contracted on this edge, i.e. the graph with the edge e deleted and the two vertices that 
it joined identihed. It is straightforward to establish that the definitions (|1.1.2|) - (p..l.3|) and 
( p..l.l|) are equivalent. The dehnition of the Tutte polynomial can be generalized from graphs 
to matroids but we shall not need this generalization here. 

An elementary consequence of the dehnition (|1 .1. 1|) , is that for a planar graph G, the 
Tutte polynomial satishes the duality relation 


T{G,x,y) = T{G*,y,x) 


where G* is the (planar) dual to G. 


(1,1,4) 
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1.2 Recursive Graphs 


Def. A recursive graph Gm is a graph that is constructed by connecting m copies of a 
subgraph H together sequentially in some prescribed manner. 


One important class of recursive graphs is provided by strips of regular lattices, such as 
the square lattice. Let us envision such a strip as extending in the longitudinal direction, 
along the x axis, and having length = m vertices, with width Ly along a transverse [y) 
axis (where no confusion should result between these axes and the x, y variables of the Tutte 
polynomial). Let Piy denote the path graph containing Ly vertices, denoted ui, ^ 2 ,..., 
Formally, one can dehne the strip of the square lattice by starting with the subgraph H = Piy 
and specifying the linkage L that connects two successive copies of H as the identity linkage, 
such that vertex Vi of the hrst P^ is connected by an edge to vertex v[ of the next Pl for i = 
1,..., Ly. This linkage will be denoted L = id. Finally, one specihes how the longitudinal ends 
are treated. Possibilities include free ends, a cyclic strip, such that one identihes the subgraph 
Hm+i with ifi, and a Mobius strip, such that one identihes Ffm+i with T(iLi), where T{H) 
means a twist, i.e., a reversal of transverse orientation. We shall denote the cyclic and Mobius 
square-lattice strips as sq{Ly, L^., FBCy, PBC^) and sq{Ly, FBCy,TPBCx), where, in 
physics nomenclature, FBC and PBC refer to free and periodic boundary conditions. One 
may also consider H = where denotes the circuit graph with n vertices. This 

is equivalent to periodic transverse boundary conditions, PBCy, in physics terminology. 
Then the square-lattice strip graphs sq{Ly, L^, PBCy, FBCx), sq{Ly, L^, PBCy, PBCx), and 
sq{Ly, Lx, PBCy,TPBCx) can be embedded in cylindrical, torus, and Klein bottle surfaces, 
respectively. For hxed Ly, we shall refer to the set of strip graphs sq{Ly, Lx, PBCy, PBCy) 
for variable Lx as a family, and so forth for other recursive families of graphs. An early study 
of chromatic and Tutte polynomials for recursive graphs is ||^. Some calculations of Tutte 
polynomials for lattice strip graphs are in ||15||-[p^ 


Another important class of recursive graphs is comprised of clan graphs |^. First, recall 
two auxiliary dehnitions: A complete graph Kr is a graph containing r vertices with the 
property that each vertex is connected by edges to every other vertex. (Clearly, \E\ = (0 
for Kj..) The join of two graphs Hi and H 2 , denoted Hi + H 2 , is the graph formed by 
connecting each vertex of Hi to all of the vertices of H 2 with edges. Then 


Def. A clan graph is a recursive graph composed of a set of m complete graphs 
Kr 2 ,..., such that the linkage between two adjacent pairs, say, Kr^ and is a join. 

Def. A homogeneous clan graph is a clan graph with the property that all of the Ky. 
are the same, say Ky. 


Def. A cyclic clan graph of length m is a clan graph with the Kyi’s arranged around 
a circle, i.e., with the identihcations Ky^^^ = Ky. (i.e., with periodic longitudinal boundary 
conditions, in physics terminology). 

Def. A homogeneous cyclic clan graph is thus a cyclic clan graph for which r* = 
r V i; we shall denote this as G[{Ky)yn, jn], where jn indicates the join linkage between each 
successive pair of Ky’s. We shall omit the qualiher “homogeneous” where it is obvious from 
the context. 
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Previously, we calculated the Tutte polynomial for the family of cyclic clan graphs of 
length m composed of K 2 ^s, i.e., G[{K 2 )m, jn] [2^. This family is equivalent to the fam¬ 


ily of cyclic strips of the square lattice of width Ly = 2 and length = m with next- 


nearest-neighbor spin-spin interactions. The asymptotic accumulation set of the zeros of the 
corresponding Potts model partition function was determined in the limit m 
We remark on some basic properties of the graph G[{Kr)m, j"^]- This has 


oo. 


\V\ = mr , |£^| = 


mr{3r — 1) 


( 1 . 2 . 1 ) 


and is a A-regular graph with uniform vertex degree A = 3r — 1 (here a A-regular graph 
is one in which all vertices have the same degree A). It is straightforward to show that 
if one cuts the cyclic clan graph G[(Kr)m, j'>^] cit one of linkages between successive K^s, 
say and and then reconnects the vertices of these two complete graphs with each 

other with a twist, one obtains precisely the same graph. Recall that a clique of a graph 
G is dehned as the maximal complete subgraph of G. For m sufficiently large as to avoid 
degenerate special cases, the clique of G[{Kr)m, jn] is K 2 r- 


1.3 Equivalence Between Tutte Polynomial and Potts Model Par¬ 
tition Function 


Besides its interest for graph theory, the Tutte polynomial has a close connection with 
statistical physics, since it is, up to a prefactor, equal to the partition function for a certain 
model of phase transitions and cooperative phenomena known as the g-state Potts model 

0,11, E 


We review this connection here since we shall use it below. 


Def. For a statistical mechanical system in thermal equilibrium at temperature T, the 
partition function (= sum over all spin conhgurations {a} ) for the g-state Potts model on 
the graph G = (V, E) is given by 


Z{G,q,v) = J2e-^ ( 1 . 3 . 1 ) 

W} 

with the Hamiltonian describing the spin-spin interaction 

H = ( 1 . 3 . 2 ) 

(u> 

where Uj = 1, ...,g are the effective spin variables on each vertex i G R; J is the spin-spin 
coupling; (3 = {kBT)~^] ks is the Boltzmann constant; and (ij) denotes the edge in E joining 
the vertices i and j in V. 

We use the notation 

K = j3J , a = , v = a — 1 (1.3.3) 

(where v should not be confused with V) so that the physical ranges are (i) J > 0 and 
hence a > 1, i.e., u > 0 corresponding to cx) > T > 0 for the Potts ferromagnet, and (ii) 
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J < 0 and hence 0 < a < 1, i.e., —1 < n < 0, corresponding to 0 < T < cx) for the Potts 
antiferromagnet. An important function in physics is the (reduced) free energy per site /: 

Def. The (reduced) free energy of the g-state Potts model on the graph G = (id, E) is 

f{{G},q,v)= lim ln[Z(G,g,u)^/l^l] (1.3.4) 

|V|^oo 


where we use the symbol {G} to denote the formal limit lim|y|^oo G for a given family of 
graphs. As discussed in [^, for physically relevant values of q and v, Z{G,q,v) is positive 
and hence it is clear which of the l/|ld|’th roots (of which there are |ld|) one should choose 
in evaluating (|1.3.4|) . However, in other regions of the space of (g, v) variables, Z{G, g, v) can 
be negative or complex, so that there is no canonical choice of which root to take in ( |1.3.4| ) 
and only the quantity | exp(/({G}, g, v))! can be obtained unambiguously. 

To show the equivalence of the Potts model partition function and Tutte polynomial, we 
shall recall the elementary theorem 

m 


Theorem 

expressed as 


The Potts model partition function on the graph G = (V, E) can be 
Z{G,q,v) = (1.3.5) 


G'CG 


where G' is a spanning subgraph of G. 

Proof This theorem is proved by observing that 

Z{G,q,v) = J2ll^^ + vSo 

G} (o) 


(1.3.6) 


and carrying out the indicated product over edges and sum over spin conhgurations. □ 

The formula ( p..3.6|) allows one to generalize g from Z_|_ to M; more generally, it allows 
one to generalize g and v to C. 

One then has the well-known and important result 


Corollary For a graph G = {V,E) the Potts model partition function Z{G,q,v) is 
related to the Tutte polynomial T{G,q,v) according to 



Z{G, g, v) = {x- - 1)I^It(G, x, y) 

(1.3.7) 

where 

x=l+^- 

(1.3.8) 


V 

and 

y = a = V + 1 

(1.3.9) 

so that 

1 

1 

(1.3.10) 


Proof This follows immediately from (|1.3.5|) and the dehnition (|1 .1. 1|) . □ 
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1.4 Chromatic, Flow, and Reliability Polynomials 

An important property of the Tutte polynomial is the fact that it is a Tutte-Grothendieck 
(TG) invariant 


Thus, consider a function / (not to be confused with the reduced 


free energy / in (|1.3.4|) ) that maps a graph G to the elements of some held K. Let fb = 
f (bridge) and fi = f(loop). It will suffice to consider a connected graph G; for disconnected 
graphs, one requires that f(Gi U G 2 ) = f(Gi)f(G 2 ). 

Def. The function / : G —y K is a Tutte-Grothendieck invariant if (1) if e G -E is a 
bridge, then f(G) = fbf(G/e), (2) if e G -E is a loop, then f(G) = fif(G/e), (3) if e G -E is 
neither a bridge nor a loop, then 


Next, we recall the 


Theorem [11 


f(G) = af(G-e) + bf(G/e) a,b^0 


f(G) = —) . 

0 a 


(1.4.1) 


(1.4.2) 


For the proof, see [111. As a consequence of this, many important polynomial functions 


of graphs that are expressible in terms of TG invariants are particular cases of the Tutte 
polynomial. To render our discussion self-contained, we review the relevant dehnitions and 
relations here. 

The hrst special case is the chromatic polynomial P{G,q), which counts the number of 
ways of coloring the vertices of G subject to the constraint that no adjacent pairs of vertices 


have the same color 30 - 34 . Then 


P(G,q) = (-,)‘<°'(-l)''''r(G,l-«,0) . 


(1.4.3) 


This result follows from (|1.4.2|) by observing that q~^P is a TG invariant with q~^P(bridge) = 
q — I, q~^P{loop) = 0, satisfying ( [1.4.1|) with a = 1, b = —1. A natural way for a physicist 
to prove this is to observe hrst that the chromatic polynomial is identical to the n = — 1 
special case of the Potts model partition function: 


P(G,q) = Z(G,q,v = -l) 


(1.4.4) 


since for this value (corresponding to the zero-temperature Potts antiferromagnet) no two 
adjacent spins can have the same value. One then uses (|1.3.7|) to obtain ( |1.4.3|) . 

Ghromatic polynomials have been of interest in mathematics for many years since Ref. 
H , owing in particular to their connection with the proper face-coloring of bridgeless planar 
graphs and early efforts to prove the four-color theorem, that a proper face-coloring of a 
bridgeless planar graph G can be accomplished using four colors. This theorem is equivalent 
to the statement that if G is a loopless planar graph, then P(G, 4) is a positive integer, 
i.e., there is a proper vertex coloring of G with g = 4 colors. Ghromatic polynomials are of 
interest for statistical physics because they are equal to the partition function of the g-state 
Potts antiferromagnet at zero-temperature (i.e., v = —1), as indicated in (|1.4.4|) above. The 
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degeneracy of spin states per site for this g-state Potts antiferromagnet on the |id| —> cx) 
limit of the graph G (usually a regular lattice with some specihed boundary conditions) is 
defined as 

W{{G},q)= lim . (1.4.5) 

|y|— >oo 

As discussed in |^, for the recursive graphs of interest here, for sufficiently large real q, 


P{G,q) is positive, so that one has a natural choice for which of the 1/\V\ roots to pick in 
evaluating ( 1.4.5|) . In terms of the singular locus B to be defined below, one defines the region 
Pi to be the maximal region in the q plane to which one can analytically continue W ({G}, q) 
from the interval q > gc({G}), where gc({G}) denotes the maximal point where B intersects 
the real axis |^]. For real q > qc, and hence for all of region Pi, there is no ambiguity 


in choosing (|1.4.5|) . However, for regions Pj not analytically connected with Pi, only the 
magnitude |hF({G},g)| can be determined unambiguously since there is no canonical choice 


of which of the 1/|H| roots to pick in evaluating ( |1.4.5|) . The ground state entropy is dehned 

5o({G},g) = fcslnlT({G},g) (1.4.6) 


as 


where /cs is the Boltzmann constant. The importance of this is that the Potts antiferro¬ 
magnet exhibits nonzero ground state entropy for sufficiently large g on a given {G}, an 
exception to the third law of thermodynamics 1^. A physical system that exhibits nonzero 


residual entropy at low temperatures is ice [^. Previous studies of chromatic polynomials 

a-m, i 


include 


We next discuss the flow polynomial and start by recalling its definition. A recent review 


IS 


Def. Consider a connected graph G = (V, E) with an orientation of each edge e G P, 
and assign to each edge a nonzero number in the additive group Zg/{0}, i.e., 1, 2,..., g — 1. 
A (nowhere-zero) flow on G is an assignment of this type satisfying the property that the 
ingoing and outgoing flows are equal mod g at each vertex. The number of such flows is 
given by the flow polynomial P(G, g). 

Remarks: As is implicit in the notation, P(G, g) does not depend on the orientation 
chosen for the edges. Moreover, one can generalize the notion of a flow to involve the 
assignment of an element of a nonzero element of any abelian group of order g, not just Zg, 
to the edges of G, but we shall not need this generalization here. For an arbitrary bridgeless 
graph it has been proved by Seymour that there is always a (nowhere-zero) 6-flow |^|, and 
it has been conjectured by Tutte that there is always a (nowhere-zero) 5-flow [|^. 

The flow polynomial P(G, g) has the following properties. It obviously vanishes on a 
bridge, 

F{bridge,q) = 0 (1-4-7) 


and satishes 


F{loop, g) = g — 1 

F(G, q) = (q — 1)P(G — e, g) if e G P is a loop 


(1.4.8) 

(1.4.9) 
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Further, 


F{G, q) = F{G/e, q) — F{G — e,q) if e G -E is not a loop . 


(1.4.10) 


Hence, from (|1.4.2|) the flow polynomial for a (connected) graph G is given in terms of the 
Tutte polynomial by 

F{G, q) = (-l)l^l-l^l+iT(G, 0,1 - g) . (1.4.11) 

An elementary result is that if G is planar, then P(G, q) = qF{G*, q), where G* again denotes 
the planar dual of G. 

Third, for a connected graph G = {V,E), consider a related graph if obtained from 
G by going through the full edge set of G and, for each edge, randomly retaining it with 
probability p (thus deleting it with probability 1 — p), where 0 < p < 1. The (all-terminal) 
reliability polynomial R{G,p) gives the resultant probability that any two vertices in G are 
connected, i.e., that for any two such vertices, there is a path between them consisting of a 
sequence of connected edges of G H- The reliability polynomial of the connected graph G 
is given in terms of the Tutte polynomial by 


P(G,p) = pl'"l-'(l - p)l^l-l'"l+'T(G, 1,1/(1 - p)) . 


(1.4.12) 


A basic property is that for 0 < p < 1, it follows that 0 < R{G,p) < 1. Note that since, 
obviously, R{G, 1) = 1 for the connected graph G, it follows that the prefactor (1—p)l-®l“l^l+i 
in (|1.4.12|) is always cancelled by an inverse factor (1 — p'j~\E\+\v\-i T{G, 1,1/(1 — p)) 
so that R{G,p) has no (1 — p) factor. Studies of roots of the reliability polynomial include 

HI, 11,0. 

Special cases of Tutte polynomials also yield, up to prefactors, other important functions, 
such as the Jones polynomial Viit) for alternating knots and links, obtained by setting x = t 
and y = -1/t Q (see, e.g. [0, |1). 


1.5 Special Valuations of Tutte Polynomials 


For certain special values of the arguments x and y, the Tutte polynomial T{G,x,y) yields 
various quantities of basic graph-theoretic interest i-0- We recall that a tree is a con¬ 
nected graph with no cycles; a forest is a graph containing one or more trees; and a spanning 
tree is a spanning subgraph that is a tree. We also recall that the graphs G that we consider 
are connected. A valuation that will be of interest below is the number of spanning trees of 
a graph G, denoted by Nst{G). From the dehnition (|1.1.1D , it is immediately evident that 
this is given by 

Nst{G) = T{G,1,1) . (1.5.1) 


A second valuation of the Tutte polynomial that will be of interest here yields the number 
of acyclic orientations. Thus, consider a connected graph G and dehne an orientation of G 
by assigning a direction to each edge e ^ E. Clearly, there are 2^^^ of these orientations. 


Def. For a connected graph G, an acyclic orientation is dehned as an orientation that 
does not contain any directed cycles. (Here a directed cycle is a cycle in which, as one 
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travels along the cycle, all of the oriented edges have the same direction.) The number of 
such acyclic orientations is denoted a{G). 

A basic result, due to Stanley, is 

a(G) = q = -l) . (1.5.2) 

Equivalently, in terms of the Tutte polynomial, 

a{G) =T{G,x = 2,y = 0) . (1.5.3) 

Other valuations include the number of spanning forests, Nsf{G) = T{G, 2,1), the number 
of connected spanning subgraphs, Ncssg{G) = T(G, 1,2), and the number of spanning 
subgraphs, Nssg{G) = T(G, 2,2) = 2l®l, where the latter relation holds since T{G,2,2) is 
the number of spanning subgraphs of G, and in counting these, one has the two-fold choice 
for each edge of whether to include it in G' or to exclude it from G'. 


1.6 Form of Tutte Polynomials for Recursive Families of Graphs 


In [|^, the following theorem was proved: 

Theorem 1 Let Gm be a recursive graph. Then the Tutte polynomial and Potts 

model partition function have the general forms 


^T,G,X 

T{Gm,X,y) = ^ CT,G,j{^T,G,jy 
i=l 

^z,G,\ 

Z(Gm,q,v)= ^ Cz,Gj{Xz,G,jy 


( 1 . 6 . 1 ) 


( 1 . 6 . 2 ) 


where 


Nt,g,x = Nz,g,\ ■ (1.6.3) 

The coefficients CT,G,j and terms \T,G,j in eq. (|1.6.1|) and the coefficients cz,G,j and terms 
Xz,G,j in eq. ( [1.6.2|) do not depend on the length m of the recursive graph. The coefficients 
cz,G,j only depend on q, and satisfy 


cr,G,j = 


CT,G,i 


X 


(1.6.4) 


CT,G,j - Cz,G,j 


(1.6.5) 


with 

See [M for further discussion. 

Remark: A related property, shown in [0, is that the Tutte polynomials T{Gm,x,y) for 
recursive graphs of different lengths m satisfy a hnite recursion relation. The special case of 
















( p..6.1|) for the chromatic polynomial was noted earlier in [^, namely that for a recursive 
graph 

P{Gm,(l)= ^ Cp,G,j{^P,G,j)^ ■ ( 1 . 6 . 6 ) 

i=i 

We next dehne the sums of coefficients cz,G,j- 

Def. Let Gm, be a recursive graph so that the structural theorems above hold. Then 

^Z,G,\ 

Gz{G) = ^ cz,G,j 

i=i 

^T,G,X 

= CT,G,j 

i=i 


= C't(G) 


(1.6.7) 


and 




( 1 . 6 . 8 ) 


^P,G,\ 

Gp{G) = ^ cp,G, 

i=i 

It was shown earlier how these sums are related to colorings of the repeated subgraph ^ 


Before proceeding with our new results, we recall the standard notation from combina¬ 
torics for the falling factorial q^r) and rising factorial 


r—1 

q(r) = 

r—1 

gh) = JJ(g + s) . 

s=0 

These satisfy the relations 

gW = (g + r _ 



rp) = rp_i) = r! . 


(1.6.9) 

( 1 . 6 . 10 ) 

( 1 . 6 . 11 ) 

( 1 . 6 . 12 ) 

(1.6.13) 

(1.6.14) 

(1.6.15) 
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1.7 Continuous Accumulation Sets of Zeros 


In the context of cyclic clan graphs, we shall sometimes use the symbol G[{Kj.)oo, j'n] for 
the formal limit of the family as m —> cx). Since T{G,x,y) is equal, up to a prefactor, to 
Z{G,q,v), we shall usually, with no loss of generality, consider the zeros in the space 
spanned by the Pott model variables {q, v) rather than the zeros in the space spanned by 
the Tutte variables {x, y). 


Def. Consider a recursive graph Gm and the limit m ^ oo. The continuous accumu¬ 
lation set of zeros of Z{Gm,q,v) as m —>• oo is denoted B. The slices of this locus in the q 
plane for fixed v and in the v plane for fixed q will be denoted Bq and B„, respectively. 

Remarks: In cases where the locus B is nontrivial, it forms by the merging together of 
the zeros as m —> cxo. This continous locus may be empty if the zeros accumulate at one or 
more discrete points. A simple example of the latter is provided by the special case of the 
chromatic polynomial of the tree graph, Z(Tm, q, —1) = P(Tm, q) = q{q — 1)™'“^. Evidently, 
this polynomial has only the discrete zeros q = 0 with multiplicity one and q = 1 with 
multiplicity m — 1, so that there is no continuous accumulation set; Bq = ^. A simple 
example of a nontrivial locus is provided by the Potts model partition function of the circuit 
graph. An elementary calculation yields Z{Gm, q, v) = {q + + (<? — l)n'", or equivalently, 

T{Cm, X, y) = y + J2^Si xK The locus B is given by the solution to the equation |g-|-n| = |n|. 
For fixed n, the locus Bq is the circle centered at g = —v with radius |n|. For fixed real q, 
the locus By is the vertical line Re(n) = —g/2. For the corresponding Tutte polynomial, B 
is the solution of the equation |a:| = 1, so that B^, in an obvious notation, is the unit circle 
in the x plane and there is only a discrete zero in the g-plane, so that By = 0. 


Def. Consider a recursive graph G^, so that Z{Gm, q, v) has the form (|1.6.2|) , and let 
q and v be given. Among the terms \z,G,j in ( |1.6.2|) the one with the maximal magnitude 
\^z,Gq \ at this point (g, v) will be denoted as “dominant”, labelled \z,G,dom.- This corresponds 
to the dominant term \T,G,j in (|1.6.1|) at the corresponding point {x,y). 

Remark: Since the Xz,G,j are distinct, for a generic (g,n), only one term Xz,G,j will 
be dominant. As m —*• oo, the reduced free energy will thus be determined only by this 
dominant term: 

/ = \^[{Xz,G,dom)^^^] (IT.l) 


where 

t= lim . (1.7.2) 

m^oo m 

As one moves to another point (g', v'), it may happen that there is a change in the dominant A, 
from Xz,G,dom. to, say, X'^ c^dom- if) indeed, this happens, then there is a resultant nonanalytic 
change in the free energy / as it switches from being determined by Xz,G,dom. to being 
determined by X'zGdom..- Hence 

Theorem 2 Consider a recursive family of graphs Gm and let m — > cx). Then B is 
determined as the solution to the equation 


|A 


Z,G,dom. 


— \^z, 


G,dom. 


(1.7.3) 
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A corollary is that although /({G}, q, v) is nonanalytic across B, the quantity | exp(/({G}, q,v))\ 
is continuous across this locus. 

Thus, for the Potts model partition function on a recursive family of graphs Gm, in 
the m —>■ oo limit, the resultant locus B forms a hypersurface in the space spanned by 
{q,v), and the reduced free energy is nonanalytic on this surface. One can determine the 
corresponding slices Bq in the q plane for hxed v and in the v (or equivalently, a) plane for 
hxed q. In our earlier calculations of Tutte polynomials and Potts model partition functions 
for recursive families of graphs this program was carried out 

Similarly, for the one-variable specializations of the Tutte polynomial, we denote the 
corresponding continuous accumulation set of zeros in the respective variables as Bq for the 
m ^ oo limits of the chromatic and flow polynomials P{Gm,q) and F{Gm,q), and Bp for 
the m ^ oo limit of the reliability polynomial R{Gm,p)- For a given {G}, as one crosses the 
locus Bq, there is a nonanalytic change in the form of the function W ({G}, q) associated with 
a switch in the dominant term Xp,G,dom.- Similar remarks hold for the flow and reliability 
polynomials. 

As a consequence of the fact that the Tutte polynomial has real coefficients and of the 
relations connecting the chromatic, flow, and reliability polynomials with the Tutte poly¬ 
nomial, it follows immediately that the set of zeros of each of these polynomials for any m 
in the respective complex planes is invariant under complex conjugation, and, in the limit 
m ^ oo, the respective continuous accumulation sets B are also invariant under complex 
conjugation. 

Since for g G and for sufficiently large v, Z{Gm,q,v) grows exponentially rapidly as 
m —>■ oo, as can be seen from the structural theorem (|1.6.2| ), the limit in eq. (|1.3.4|) dehnes 
a hnite reduced free energy /({G}, q, v) as |Id| —> oo. Similarly, for sufficiently large positive 
q, P{Gm, q) grows exponentially as m —>■ oo, and the dehnition in eq. (|1.4.5| ) yields a hnite 
W{{G},q). 

Similar remarks hold for the how polynomial F{Gm, q) and reliability polynomial, and 
we are thus led to introduce the dehnitions 

Def. Let Gm be a recursive family of graphs and consider the limit m ^ oo. For 
sufficiently large real q, F{Gm, q) is real and positive, and we dehne 

(/.({G},g)= lim F(G™,g)i/l^l (1.7.4) 

y|^oo 

with the canonical choice of l/|I7|th root so that 0({G}, q) is real and positive. The function 
0({G},g) measures the number of (nowhere-zero) g-hows per vertex in this limit. 

Def. Let Gm be a recursive family of graphs and consider the limit m ^ oo. For 
0 < p < 1, we dehne 

p{{G},p) = lim R{Gm,pf^''^'' ■ (1.7.5) 

|V|—>oo 

The function p({G},p) is a measure of the connectivity, normalized to be per vertex, in this 

limit. 


27 
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1.8 Chromatic Polynomial for Cyclic Clan Graphs 

Since this paper includes new results on the Tutte polynomial of cyclic clan graphs, it is useful 
to recall that the special case of the chromatic polynomial for the clan graph G[{Kr)m, jn] 
was calculated by Read |@; in our present notation, it is 


where 


^P,r,d — ( — 




Q 

d-1 


Note that 


and 


m.)jdi],g) ^ ^ /^d(-^P,r,d) 

(1.8.1) 

d=0 


^ =( !)%)(? r c?)(r-d) 

Q(r+d) 

(1.8.2) 

in g of degree d given by 


/iO = 1 

(1.8.3) 

(,-2d+l) 

d\ 

(1.8.4) 

^P,r,0 (g ^)(f) 

(1.8.5) 

Ap,r,r = (-l)V! 

(1.8.6) 

= (—l)’'“^r!(g — 2r + 1) . 

(1.8.7) 


The first few coefficients Hd beyond po are pi = g — 1, /i 2 = (l/2)g(g — 3), ps = (l/3!)g(g — 
l)(g —5), and = (l/4!)g(g —l)(g —2)(g —7). For our discussion below it will be convenient 
to extract the factor of q that is present in for d >2 and define 


fid = g Vd for d>2 


The sum of the coefficients is 


Cp{G[{Kr)m,jn]) = 


( 1 . 8 . 8 ) 


(1.8.9) 


Parenthetically, we note that Read actually proved a more general result, namely the chro¬ 
matic polynomial for an inhomogeneous cyclic clan graph where the ri, r 2 ,..., are, in 
general, different [^|. This chromatic polynomial does not have the form (|1.6.6 ) since the 
general inhomogeneous cyclic clan graph is not a recursive graph. Biggs and coworkers have 
also obtained a number of results on chromatic polynomials of recursive families of graphs 
comprised of m copies of for arbitrary linkage L, G[{Kr)m, L] |^, ^ |8^. These are clan 
graphs if and only if the linkage is a join, L = jn. 
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1.9 Organization 

Having completed the discussion of relevant definitions and background, we now proceed 
to present our new results. In next section, we shall determine the structure of the Tutte 
polynomial for the cyclic clan graph jhi]. These general results will be illustrated by 

our previous calculation of the Tutte polynomial for the case r = 2, and here we shall present 
a calculation for the case r = 3 (for arbitrary m). For the consideration of the limit m —> cx), 
it is convenient to work with the equivalent Potts model partition function Z{G[{K^)mi jn], 
and we shall give some results for the B for this function. Next, we present calculations of 
two special cases of Tutte polynomials, namely, flow and reliability polynomials, for cyclic 
clan graphs and discuss the respective continuous accumulation sets of their zeros in the 
limit m ^ oo. We then prove two theorems that determine the number of spanning trees 
on G[{Kr)m) jn] and G[{Kr)mGd]- Finally, we report calculations of the number of acyclic 
orientations for strips of the square lattice and use these to obtain an improved lower bound 
on the exponential growth rate of the number of these acyclic orientations. 


2 Structural Relations for Tutte Polynomials of G[{Kr)rm jn 

In [^, relations were proved connecting the coefficients cr,Gj = cz,G,j in the Tutte poly¬ 
nomial and Potts model partition function for a recursive graph G with the corresponding 
coefficients in the chromatic polynomial P{G,q). For sufficiently large integral q these coef- 
hcients are determined by the multiplicities of the the terms \T,G,j or equivalently Xz,G,j as 
eigenvalues of an appropriately dehned coloring matrix (for the chromatic polynomial case, 
see [^). From this correspondence, it follows that the types of coefficients that enter in the 
chromatic polynomial are the same as those that enter in the Tutte polynomial. Then 

Lemma 1 

1 

a; — 1 


nT{r,d) 


T{G[{Kr)mGn\,x,y) = TtE 

d=0 j=l 


( 2 . 1 ) 


Equivalently, 


r nT{r,d) 

Z[G[{Kr)mi j'^]^ Qi'^) “ h'ci ^ ^ {,Xz,r,d,j) 

d=0 i=l 


( 2 . 2 ) 


Proof This follows from eq. (|1.6.5|) together with (|1.8.1|) . □. 
Using (|1. 3.101) , (|1.8.4|) , and (|1.8.8| ), one can re-express (^Tj) as 


T{G[{Kr)mGn],x,y) = 


X 


nj'ir^d) r nj'(r^d) 

E E + ( 2 /~ 1) EF'i E (2-3) 


d=0,l 


i=i 


d=2 


i=i 


where the second sum is understood to be absent for r < 2. 
We next observe that 
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Lemma 2 


(2.4) 


Cz(G|(/4)™,jn]) =CT(G|(AV)„,jn])) = ^ . 

Proof Because the linkage between each repeated Kr is a join, one can freely permute 
the labels on the vertices of each K^. From the sum over states in q, v) together with the 
relation (|1.3.7|) , it follows that the total multiplicity of colorings is enumerated by assigning 
colors freely to each vertex of a given Kr, taking into account the above equivalences. But 
this number is precisely g^^^/r!. For example, for the case r = 3 and q = 3, one has the 
colorings (111), (222), (333), (112), (113), (221), (223), (331), (332), and (123). We note the 
equivalences (112) = (121) = (211) and so for for colorings with two repeated colors, and 
(123) = h(123) where h is an element of the permutation group S'3. Hence, the total number 
of inequivalent colorings is 10 = 3®/3!. □ 

We now present our structure theorem; 

Theorem 3 For G[{Kr)m-i j^] th® numbers nT{r,d) = nz{r,d), 0 < d < r, are deter¬ 
mined by 

nrir, r) = 1 (2.5) 

ny(r,0) =nr(r, 1) = 2’'-^ (2.6) 

together with the recursion relation 


nT^r + l,d) = riTir, d) + riTir, d — 1) for 2 < d < r . (2.7) 


For d > r, riTir, d) = 0. 

Proof From (|2.4|) we have 

CT{G[{Kr)mdn])) = '^nT{r,d)fid = ^ • (2.8) 

d=0 


We differentiate this equation r times to get r -|- 1 linear equations for the r -|- 1 unknowns 
?7,(r(r, d), d = 0,1,..., r. Solving this equation yields the results in the theorem. □ 

Corollary 1 For G[(iFr.)m, jh]. 


nrir, d) 



(2.9) 


Proof These results follow in a straightforward manner by explicit solution of the set 
of r -|- 1 linear equations used in the proof of Theorem 3. □ 

For example, one has the special cases 

nT{r,r — 1) = r (2-10) 
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nT(r, 2) = 2'’ ^ — 1 for r > 2 
nT{r, 3) = — r for r > 3 


and so forth. 

In Table 1 we show the nnmbers nT{r, d) and Nrp^r,x for 1 < r < 8. 
Dehning iV'r j,,a as 

r 

NT,r,x = '^ririr, d) 
d=0 


we hnd 

Theorem 4 

A'7-,,A = (r+3)2’-V 


( 2 . 11 ) 

( 2 . 12 ) 


(2.13) 


(2.14) 


Proof From (p.6|) and the recursion relation Q for the set of coefficients nT{r, d) 
with a given r (the r’th row of Table 1), it follows that 


r 

Y, Mr. d) = 2NT,r-i,x - 2"-' . (2.15) 

d=i 

Adding the hnal term nT{r, 0) and using (|2.6|) , we have 

M,r,X = 2NT,r-l,X + 2'’ ^ . (2-16) 


Combining this recursion relation for Nrp^r,x with the initial value = 2 yields the 

theorem (|2.14|) . □ 

Remark: From our explicit calculation of T{G[{K 3 )m,jn\,x,y), we have found that, in 
contrast to the situation for cyclic strips and related self-dual planar strip graphs 
for a given r, some XT,r,d with different d may be equal to each other. This does not happen 
for r = 1 or r = 2; for r = 3, we hnd (see eq. ( |2.25|) below) that At, 3 , 0,1 = At, 3 , 2 , 3 - Hence the 
total number of distinct XT,r,d,/s for r = 3 is 11. 

Our result (p.5|) shows that for d = r, there is a unique term XT,r,d=r- We have 

Corollary 2 For G[{Kr)m, jn], 


XT,r,r = r\ . 


(2.17) 


Proof Given the uniqueness, nT{r,r) = 1, this follows from the analogous term in the 
chromatic polynomial, eq. ( p..8.6|) in conjunction with the structural equation (|2.1| ) and the 
general relation (|1.4.3|) . □ 

Since XT,r,r is unique, we omit the hnal index j, setting XT,r,r,i = XT,r,r- 
The previously known calculations of T{G[{Kr)m,jn\,x,y) illustrate our general struc¬ 
tural theorems. For r = 1, the family G[{Ki)m, jn] is identical to the circuit graph with m 
vertices, Gm, and an elementary calculation gives 


m—1 


T(G^,x,y) = y+Y 
i=i 


x^ = 


X 


+ hi] — 


X 


[x™' + xy — X — y] 


(2.18) 
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Table 1: Table of numbers nT{r,d) and their sums, NT,r,\ for G[{Kr)m,j'n]. Blank entries are zero. 


r 1 d ^ 

0 

1 

2 

3 

4 

5 

6 

7 

8 

R^T,r,\ 

1 

1 

1 








2 

2 

2 

2 

1 







5 

3 

4 

4 

3 

1 






12 

4 

8 

8 

7 

4 

1 





28 

5 

16 

16 

15 

11 

5 

1 




64 

6 

32 

32 

31 

26 

16 

6 

1 



144 

7 

64 

64 

63 

57 

42 

22 

7 

1 


320 

8 

128 

128 

127 

120 

99 

64 

29 

8 

1 

704 


SO that At, 1,0 = X and, in accord with (|2.17|) , Ar,i,i = 1- 

For r = 2, the Potts model partition fnnction and eqnivalent Tntte polynomial for the 
family G[{K 2 )m, jn] were calcnlated and analyzed in One has 


r(G|(Ay,„, jn], I, y) = — + y, + /12(At,2.2)’”1 (2.19) 




i=i 


where 


A 


T,2,0,j 


y' 


+ 2t/2 + 3t/ + + 3x + 2 ± ^/R. 


20 


for j = 1, 2 


R 20 = 4 + 12x + 12y + 22xy + 13a;^ + 21y‘^ + + 20t/^ + 16xt/^ 

+ lQx‘^y + + lOt/^ — Ax^y"^ — 2y^x + 4t/® — 2x^y^ + y^ 


A 


1 r 




y' 


+ 2y‘^ + 3y + 2x + 4 i R 21 


for j = 1, 2 


2R21 = 16 + 16a: + 32y + 4x'^ + 12xy + 33y^ + 20y^ - 4xy^ + lOt/^ + 4y^ + y^ (2.23) 


( 2 . 20 ) 


( 2 . 21 ) 

( 2 . 22 ) 


and At, 2,2 = 2, in accord with (|2.17|) . Thns, 77,72(2,0) = 7772(2,1) = 2 and 7772(2,2) = 1, as in 
Table |^. 

Using a systematic application of the deletion-contraction property, we have calcnlated 
the r = 3 case, i.e., T(G[(iF3)m, jn),a:,?/). Onr general strnctnre theorem yields 7772(3,0) = 
7772(3,1) = 4, 7772(3, 2) = 3, and 7772(3, 3) = 1, and we have 


T(Gl(K3)rn,jn],x,y) = 


1 


X — 1 


y^(AT,3,0,j)"' + hi y^(AT,3,l,j)' 


J=1 


i=i 
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(2.24) 


3 

+ /^2 ,)™ + /i3(AT,3,3 

i=i 


m 


where 


At,3,0,1 — At,3,2,3 — ^y'^iy + 1) 


A 


T,3,2,j 


3 r 
2 


y' 


+ 3y'^ + 2x + 6y + 8± a/^ 


(2.25) 

for j = 1, 2 (2.26) 

/?32 = 64 + 32x + 112y + 4a;2 + 2Axy + lOOy^ + Axy^ + 52y^ - Axy^ + 21y^ + 6y^ + / (2.27) 

and At, 3,3 = 3! in accord with CT - The other XT, 3 ,d/s are the solutions of quartic and 
cubic equations that are somewhat lengthy and hence are given in the appendix. 

In our previous works presenting other calculations of Tutte polynomials for recursive 


families of graphs []I^ [T^ we have given detailed plots of the continuous accumu¬ 

lation set B in the q plane for the m —> cx) limit of the Potts model partition function, for 
various values of v and B in the v (or equivalent a) plane for various values of q. Here we 
shall restrict ourselves to one interesting comparison. 

We begin by working out the locus Bg for the v = —1 case where the Potts model partition 
function reduces to the chromatic polynomial. Here, using (|1.8.1| ) and (|1.8.2|) , we hnd for 
the m —>■ cx) limit of the general G[{Kr)m, j'^] clan graph that this continuous accumulation 
set is the union of r circles 3 = 1 , •••, 


B = 




(2.28) 


where is the solution of the equation |g — r — j -f 1| = r — j -|- 1, viz., 

Crj : q = r+ j —l + {r— j + l)e*® , 0<6'<27r, l<j<r 

The j’th circle Crj intersects the real q axis at the points q = 2{j — 1) and 

q = qci{{Kr, K 2 r)}) = 2r 


(2.29) 


(2.30) 


Thus, these r circles osculate, i.e., coincide with equal (vertical) tangents, at g = gc, which 
is thus a tacnodal multiple point, in the terminology of algebraic geometry. 

The circles comprising B separate the g plane into the following r -|- 1 regions Rrj (where 
ext(Cr,j) and int(Cr,j) denote the exterior and interior of the circle Crj: 

Rr^i = ext(Cr,i) , (2-31) 


and 


Rrj = int(Cj,,j-i) Pi ext(Cr-,j) for 2 < j < r , 


R^ ^_j_ ^ int (fl^ T’) * 


(2.32) 

(2.33) 
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In region -Rr.i, the W function is 


IV = (Ap,r,oy/^ for geRr,i (2.34) 

In other regions, only the magnitude \W\ can be determined unambiguously, and we find 

\W\ = for q E Rrj , 2 < j < r + 1 . (2.35) 

We now specialize to the family G[{K^)m, jn]. For this case, if u = —1, the locus Bq 
consists of the union of the three circles : |g — 3| = 3, '■ I? ~ 4| = 2, and 

^ 3,3 • I*? ~ 5| = 1 which osculate at qc = 6. In Fig. |l| we show this locus, together with 
chromatic zeros calculated for the length m = 20. For this length, the zeros (except for the 
discrete real zeros at g = 1, g = 3, and g = 5) lie reasonably close to the asymptotic curves 
comprising the locus B. 

Using our new calculation oiT{G[{K^)mG'V\,x, y) and the equivalent Z{G[{K 3 )m, jn], g, v) 
we next show in Fig. |] the locus Bq for the illustrative value v = —0.9, together with par¬ 
tition function zeros calculated for the length m = 20. Again, except for the real zeros at 
g = 1 and g = 3, the zeros lie reasonably close to B for this value of m. Several interesting 
differences are evident: (i) increasing v from —1 to u = —0.9 removes the osculation point at 
g = 6 and replaces it with two complex-conjugate pairs of T-intersection points (such points 
were previously encountered in many loci B, e.g. [0); (ii) the right-most portion of Bq 
moves leftward, from g = 6 at u = —1 to g ~ 5 for v = —0.9; and (iii) the left-hand portion 
of B passes through g = 0 in both cases. As with W{{G},q) in the v = —1 case, different 
analytic forms for the free energy apply in the different regions bounded by the portions of 
B for general v. 


3 Flow Polynomials for G[{Kr)m^ j^] 

We hrst note the 

Corollary 3 


r nT{r,d) 

F{GHKX.,M,q) = E ET.rAi(x ^ 0,y = 1 - q)r . (3.1) 

d=0 j=l 


Proof This follows immediately from combining eqs. ( [1.4. 11 ) and (^) and using 

(EUD- ° 

It will be convenient to write (|3.1|) simply as 


r nT{r,d) 

F{G[{KrUjn],q) = E 

d=0 j=l 


(3.2) 
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Figure 1: Singular locus B in the q plane for the m ^ oo limit of the Potts model partition function 
Z{G[{K 3 )m, jn\,q,v) for v = —1, where this becomes the chromatic polynomial. For comparison, chromatic 
zeros for m = 20 (i.e., \V\ = 60) are also shown. 




Figure 2: Singular locus B in the q plane for the m ^ oo limit of the Potts model partition function 
Z{G[{K 3 )m, jn],q,v) for v = —0.9. For comparison, zeros of this partition function are shown for m = 20 
(i.e., |P| =60). 


20 



where 


Since G[{Ki)m,jn] 
tary: 


^F,r,d,j — ( — 1)^^^'’ ^^''^XT,r,d,j ■ (3.3) 

Cm, it follows that F{G[{Ki)m,jn],q) = F{Cm,q)- This is elemen- 

F{Cm,q) = q-l- (3.4) 


The flow polynomial for G[{K 2 )m, jn] can be obtained from our previous calculation of 
the Tutte polynomial for this family 


2 nT{2,d) 

F{G[{K2)m,jnlq) = E 

d=0 j=l 


where nr(2, 0) = 2, nr(2,1) = 2, n'r(2, 2) = 1, and 

1 


A 


F,2,0,j 


2 L 


(g-2)(g2-3g + 4)±/R 


F20 


for j = 1, 2 


with 


A 


Rf 20 = {q-‘2){q- 3){q^ - 5q^ + Uq^ - 20q + 12) 
1 


F,2,l,j 


q^ - 5g2 + lOg - 10 ± ^/R 


F21 


for j = 1, 2 


where 

and 


= g® - lOg^ + 45g^ - 120g^ + 208^^ - 224g + 116 . 


Af,2,2,1 = Xf^2,2 — ~2 

The flow polynomials for the hrst few values of m are 

F(G[{R\)i,jn],q)^{g-lf{g-2) 


F{G[(K 2 ) 2 ,jn],q) = (q-l)(g-2) (q - 5q +12q - 16j + 10) 


(3.5) 

(3.6) 

(3.7) 

(3.8) 

(3.9) 

(3.10) 

(3.11) 

(3.12) 


F(G|(ii'. 2 ) 3 ,iii|,(j) = (q-l)(q - 2f(q^ - 4(/ + 5)(j® - 6/ + 18(/® - Mq^ + 31q - 28) (3.13) 


It is elementary to show that in general, F{G[{K 2 )m, j'n], q) has the factor (g — l)(g — 2). 

In Fig. H we show the continuous accumulation set B of the zeros of F{G[{K 2 )m,3'^,q) 
as m —>■ cx). Some curves on B extend to complex inhnity, so that B is noncompact, in the 
g plane. Hence, it is also convenient to display B in the plane of the variable 



(3.14) 
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The locus B separates the q plane (or equivalently, the 1/q plane) into several regions, 
including four regions that contain intervals of the real axis, together with two pairs of 
complex-conjugate regions. Six curves forming three complex-conjugate pairs lying on the 
locus B extend inhnitely far from the origin of the q plane. In the u plane, these curves (taken 
in their totality) have a multiple intersection point, in the sense of algebraic geometry, at 
the origin. We introduce polar coordinates, letting 

u = . (3.15) 

We have 

Corollary 4 


On = for 0 < j < 5 (3.16) 

6 

i.e., 9u = ±7r/6, ±7r/2, and ±57r/6. 

Proof To prove this, we note hrst that the two terms that are dominant in the vicinity 
of M = 0 are 2 , 0,1 and 2 , 1 , 1 - We extract a factor of q^ from these, dehning \F, 2 ,d,j = 
Q~^^F, 2 ,d,j, express the XF, 2 ,d,j in terms of u, and carry out a Taylor series expansion of these 
in the vicinity of u = 0, to hnd 

Af, 2 ,o,i = 1 ~ 5m-|- lOu^ — -|-0 (m"^) as m —>• 0 (3-17) 

Af, 2 ,i,i = 1 — 5m-|-IOm^ — IOm^-|-0(m"^) as m —>• 0 (3.18) 

The equation of the degeneracy of magnitudes of leading terms in the vicinity of m = 0, 
|Af, 2 ,o,i| = |Af, 2 ,i,i|, yields the condition 

|m|^ cos 36*„ = 0 for |m| << 1 (3.19) 


implying that cos36'„ = 0, which yields the result (p.l6|) . □. 

The locus B crosses the real axis at three points. The largest of these we denote q^, this 
has the value 

qc{F{G[{K2)oo,jn]) = 4.079828... (3.20) 

The other two crossings occur at 

qi = 2.7760454... (3.21) 

and 

q = 2 . (3.22) 

The point qc is the larger of the two real solutions of the equation 


- 7q^ + 13q^ - g - 14 = 0 


(3.23) 


which is the equation of degeneracy in magnitudes of dominant terms |Ai?^ 2 ,o,i| = |Af, 2 ,i,i|- 
The point gi is the real solution of the equation of degeneracy in magnitudes of dominant 
terms |Af, 2 ,i,i| = \Xf, 2 , 2 \, 

g=^ - 4g2 + 7g - 10 = 0 . (3.24) 
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Figure 3: Singular locus B in the q plane for the m ^ oo limit of the flow polynomial F{G[{K 2 )m,j'n],q). 
For comparison, zeros of the flow polynomial F{G[{K 2 )m,j'n],q) are shown for m = 30 (i.e., \V\ = 60). 
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Figure 4: Singular locus B in the u = 1/q plane for the m ^ oo limit of the flow polynomial 
F{G[{K 2 )m,j'n],q). For comparison, zeros of the flow polynomial F(G[(iF 2 )m, ju], g) are shown for m = 30 
(i.e., Ifol = 60). 



The four regions that contain intervals of the real axis, together with the respective intervals, 
are: (i) i?i, q > Qc, (h) R2, Qi < q < qc, (hi) R3, 2 < q < qi, and (iv) R4, q < 2. The 
two complex-conjugate regions are denoted R^^R^ and Rq^RI] in Fig. as one moves 
counterclockwise from Ri into the upper half g-plane, one traverses i? 5 , containing the point 
q = 2 + 2 i, and then Rq, containing q = 2 i. As is evident in Fig. the density of zeros is 
rather different along different portions of B. 

In region i?i, with an appropriate choice of branch cuts, the dominant A is Ai?^ 2 ,o,i so that 

0 (G'[(iF 2 )oo,jA],g) = q e Ri . (3.25) 

In region R2, the dominant A is Air^ 2 ,i,i so that 

\(j){G[{K2)oc,jn],q)\ = |Air, 2 ,i,i|^'^^ for q e R2 ■ (3.26) 

(Recall that in regions other than Ri, only |0| can be obtained unambiguously.) In region 
R3, ^F,2,2 is dominant and 

\(j){G[{K 2 )oc,jn],q)\ = V2 for g G R 3 . (3.27) 

In region R 4 , with an appropriate choice of square root branch cut, Xf,2,i,2 is dominant and 

|0(G[(iF2)oo,jn],g)| = |Ai., 2 ,i, 2 r/' forgGR4. (3.28) 

In R 5 and Rq, Xf,2,i,2 and Ai 7 ’, 2 ,o,i dominate, respectively. We also note that there are several 
T intersection points on B. 

From our calculation of T(G[(A' 3 )m, jn],x,?/) presented here, one can obtain the corre¬ 
sponding flow polynomial F(G[(iF 3 )m, jn], g). Since the expressions are rather lengthy, we 
omit the details. 


4 Reliability Polynomials for G[{Kr)rm jn] 


We proceed to discuss reliability polynomials for G[{Kr)m, j^]- Since G[(iFi)m, jn] is the 
circuit graph Gm, the result is elementary: 

R{Gm,p)=p'^~^[m{l-p)+p] . (4.1) 


Note that when one sets x = 1, the two XT,r=i,d,j simplify: At, 1 , 0,1 = x becomes equal to 
At,i,i,i = 1, so that Hence, the total number of distinct terms AT,r=i,A = 2 in the Tutte poly¬ 
nomial is reduced to Nii^r=i,x = 1 for the reliability polynomial. This reduction generalizes 
to higher r. The reliability polynomial R{Gm,p) evidently has m — 1 zeros at p = 0 and one 
other zero, at p = m/(m — 1). As m ^ 00 , this other zero approaches 1 from above. 

The reliability polynomial for G[{K 2 )m, j'^] can be obtained from our previous calculation 
of the Tutte polynomial for this family of graphs [^]. Setting x = 1, we hnd that two of the 


hve A’s in the full Tutte polynomial become equal to two others: 

At, 2 , 0,1 = At, 2 ,i,i = - {y + 2)(p^ + 3) -|- x/R2 


(4.2) 
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A 


T, 2 , 0,2 — At, 2 , 1,2 


{y + 2 ) ( 2 /^ + 3 ) — 


where 


= 36 + AAy + 33y^ + IQy^ + + Ay^ + y^ . 


(4.3) 


(4.4) 


Hence, Nt, 2 ,\ = 5 is reduced to Nii, 2 ,x = 3 for the reliability polynomial. Setting y = 1/(1—p) 
and inserting the appropriate prefactors according to (|1.4.12| ), we obtain 


where 


R{G[{K 2 )m,jn],p) = {a 2 ,ir + K 2 )”^ - 7 t(« 2 , 3 )™ + 


mp-\l-p)U (« 2 ,i)”^-'I^ 2,1 + (« 2 , 2 )™-'/^ 2,2 


“2,j - - 


2 L( 3 - 2 p)( 4 - 6 p + 3p2)±/^] for j = l ,2 
R 2 = 36p® - 260?“^ + 793p'‘ - ISOSp® + 1236p^ - 640p + 144 


(4.5) 

(4.6) 

(4.7) 

(4.8) 

02,3 = 2(1-pf. (4.9) 

As examples, skipping the degenerate case m — 1 , the reliability polynomials for G\(K 2 ) 2 ,jn] 
for the lowest two cases m = 2, 3 are 

R{G[{K2)2,jn],p) = p^(2 -p)(6 p® - 44/ + 139p^ - 242^^ + 246^^ - 140p + 36) (4.10) 


-^2,i - 


3± 


(36 — 88 p + 69p^ — 18p^ 


x/R: 


2p 


for j = 1 , 2 


and 


R{G[{K2)3,jn],p) = p®(120p^° - 1440^^ + 7830p® - 25440p^ + 54780p® - 81840p® 


+86110p^ - 63195p3 + 31080p2 - 9300p + 1296) 


(4.11) 


In Fig. H we show the continuous accumulation set B of the zeros of R{G[{K 2 )m, jn],p) 
in the limit m ^ oo. This locus is comprised of a self-conjugate arc that crosses the real axis 
at p = 3/2 and terminates at approximately p = 1.320 ± 0.309i, together with a complex- 
conjugate pair of arcs. The upper arc extends between endpoints at 0.790 -|- 0.508i and 
1.501 -|- 0.4627 The six endpoints of the arcs occur at the zeros of i? 2 p, where there are 
square root branch point singularities in some of the A’s. This locus B does not separate 
the p plane into different regions. In passing, we observe that our results are in accord with 
the conjecture that for an arbitrary (connected) graph G, the roots of R{G,p) lie in the 
disk Ip — 1 | < 1 . 
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Figure 5: Singular locus B in the complex p plane for the m —> oo limit of the reliability polynomial 
R{G[{K2)m, jn],p). For comparison, zeros of the reliability polynomial jn],p) are shown for 

TO = 30 (i.e., |F| = 60 ). 
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5 Acyclic Orientations for Recursive Families of Graphs 

5.1 General 


We recall from the introduction the definition of acyclic orientations and the relations ( |1.5.2| ) 
and (|1.5.3|) expressing the number of acyclic orientations a{G) in terms of a valuation of the 
chromatic or Tutte polynomial. First we note that for the recursive families of graphs of 
interest here, namely cyclic or cylindrical strips of regular lattices, a{Gm) grows exponentially 
as m —> cxo. This follows from the relation ( |1.5.2|) and the structural result ( p.. 6 . 6 |) . This 
motivates the dehnition 


Def. Consider a recursive family of graphs Gm and the limit m —>■ oo. Then we define 


= .Jim . 

|y 1^00 

(5.1.1) 

The number of acyclic orientations can also be generalized |96 

to as(G), given by 

ttsiG) = (—g = —s) for s G Z 4 

(5.1.2) 

with ai(G) = a{G). 

Def. Let Gm be a recursive family of graphs, such as strips of regular lattices of length 
m, and consider the limit m —>■ 00 . Then 

= .Jim . 

If 1^00 

(5.1.3) 

Hence, 

ad{G}) = fF({G},g=-s) . 

(5.1.4) 


and, in particular, a({G}) = W{{G},q = —1). 

In our previous studies of chromatic polynomials and the resultant singular loci B for 
a variety of recursive families of graphs including strips of regular lattices, we found that 
B never intersects the negative real q axis, and, indeed, the points g = —s for s G are 
in region Ri. Hence, for the recursive families of interest here, there is no ambiguity in 
which l/|l/|’th root to choose in evaluating (|1.4.5|) at the points q = —s, s G Z+. (For 
families with noncompact B, as is evident from the plots of respective loci B presented in 
15^, |62|, 1^, Q, although B does not intersect the negative real axis, it can separate the 
point q = —s from the region i?i; we do not consider such families here.) Because the points 
g < 0 are in region Ri for the strip graphs of regular lattices, we can apply a previous 
result to infer that a^dG}) is independent of the longitudinal boundary conditions. 
Thus, for example, as({G}) is the same for the m ^ oo limit of the square-lattice strip 
graphs sq{Ly, FBGy, FBG^), sq{Ly, L^, FBGy, PBG^), and sq{Ly, FBGy,TPBG^), 


and therefore we shall omit the longitudinal boundary condition in listing the values that we 
obtain in these cases. Henceforth, we shall focus on the s = 1 case, i.e., acyclic orientations. 





















5.2 a{G) and a({G}) for Strips of the Square Lattice 

We first list some elementary results. For a tree graph, 

a(T„) = . 

For a circuit graph, 


whence 

Another easy result is 


a(a) = 2^-2 

«({^}) = = a{sq,Ly = 1) = 2 


a{sq, Ly = 2, Lj. = m, FBCy, FBCx) = 2-7™ 

From P{sq[Ly = 2, = m, FBCy, {T)PBCx], q) calculated in |0, we find 

a(sg, Ly = 2,L^ = m, FBCy, PBC^) = 7”* - 2 ■ 4”^ - 2™+^ + 5 

a{sq, Ly = 2,L^ = m, FBCy, TPBC^) = 7™ - 2 ■ 4™ + 2”^+^ - 1 
From any of ( ^.2.41) , (|5.2.5|) , or ( ^.2.61) , it follows that 

a{sq, Ly = 2, FBCy) = = 2.645751... 

From P{sq[Ly = 3,Lx = m, FBCy, PBCx], q) calculated in |^, |^, we find 

a{sq, Ly = 3,Lx = m, FBCy, PBCx) = -13 + 5(2™ + 3™ + 5™) - 2 ■ 9™ - 2^, 

27 + yMA”' /27-\/M\ 


(5.2.1) 

(5.2.2) 

(5.2.3) 

(5.2.4) 

(5.2.5) 

(5.2.6) 

(5.2.7) 




+ 


(5.2.8) 


and from P{sq[Ly = 3,Lx = m, FBCy,TPBCx], q) calculated in we find 

a{sq,Ly = 3,Lx = m,FBCy,TPBCx) = 5 - 2™ + 3™ - 5™ + 2 ■ 9™ - 2A„ 


27 + Vm\"^ 


27- \/4^\ 


with 


i=i 

where the Xsqpxysj, j = 1, 2, 3, are the roots of the equation 

+ 23^2 + 134 ^ + 202 = 0 


1 (5.2.9) 


(5.2.10) 


(5.2.11) 
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From either of (^.2.81) or (|5.2.9|) we have 


a{sq, Ly = 3, FBCy) = 


27 + vW \ 


= 2.903043... 


(5.2.12) 


we 


From the chromatic polynomial P{sq[Ly = 4, FBCy, FBC^, q) calculated in [K 
obtain 

a{sq, Ly = 4, FBCy) = {Xsq.yA,ma.f'^ = 3.040731... (5.2.13) 

where \sqxyA,max IS the largest root of the cubic equation 


- 105^^ + 1747^ - 6758 = 0 


(5.2.14) 


From chromatic polynomials calculated in and we obtain the values of a{sq, Ly, FBCy) 
for Ly up to 8. The results are listed in Table ||. 

We next consider strips of the square lattice with cylindrical boundary conditions. An 
elementary calculation yields 


a{sq, Ly = 3,Lx = m, FBCy, FBCx) = 6 ■ (34)' 


(5.2.15) 


whence 

a{sq, Ly = 3, FBCy) = ( 34 )^/^ ^ 3 . 2396 II 8 ... (5.2.16) 

From the chromatic polynomial P{sq[Ly = A, L^, FBCy, FBCx], q) calculated in we 
obtain 

a{sq, Ly = 4, FBCy) = I - - - 1 

From the chromatic polynomials P{sq[Ly, Lx, FBCy, FBCx], q) for Ly = 5,6 calculated in 
^ we compute 


= 3.3944509... 


(5.2.17) 


a{sq, Ly = 5, FBCy) = 


/527+ V200585 \ 

V 2 


= 3.4481257... 


and 


(5.2.18) 

(5.2.19) 


a{sq,Ly = 6, FBCy) = = 3.4705457... 

where Xsqpxpy6,max IS the largest root of the equation 

- 2049^^ + 547805^^ - 36633324^^ ^ 639262524^ - 2756653440 = 0 . (5.2.20) 

We have also computed a{sq, Ly, FBCy) for Ly up to 12 using chromatic polynomials cal¬ 
culated in . The results are listed in Table 

A basic property of a{sq, Ly, BCy), where BCy denotes the transverse boundary condi¬ 
tion, (free, FBCy or periodic, FBCy) is that, owing to the fact that one can calculate this 
quantity (and indeed, more generally, the chromatic or Tutte polynomial from which it is 
derived by a special valuation) by transfer matrix methods, it follows that a{sq, Ly, BCy) 
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Table 2: Values of adG}) for the infinite-length limits of strip graphs of the square lattice with width Ly 
vertices and free (F) or periodic (P) transverse boundary conditions. 


Ly 

BCy 

q;({G}) 

1 

F 

2 

2 

F 

2.646 

3 

F 

2.903 

4 

F 

3.041 

5 

F 

3.126 

6 

F 

3.185 

7 

F 

3.227 

8 

F 

3.259 

3 

P 

3.240 

4 

P 

3.394 

5 

P 

3.448 

6 

P 

3.471 

7 

P 

3.481 

8 

P 

3.487 

9 

P 

3.490 

10 

P 

3.491 

11 

P 

3.492 

12 

P 

3.493 
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is a monotonically increasing function of Ly A similar monotonicity holds for 

a{tri, Ly, BCy) for strips of the triangular lattice (see further below). We also observe that 
for a given strip width Ly, the value of a for periodic transverse boundary conditions is 
greater than the value for free transverse boundary conditions. This is easily understood 
since the transverse boundary effects are reduced when one uses periodic transverse bound¬ 
ary conditions, and hence the resultant a should be closer to that for the inhnite square 
lattice. Let us denote = lim^^^oo SCy). Given the monotonically increasing 

property of a{sq, Ly, BCy) as a function of Ly, our results show that agq > 3.49287 (rounded 
off to 3.493 in the table) for this lattice. This value is a slight improvement of the recent 


lower bounds agq > 22/7 = 3.1428... and agq > 3.41358 |^. Since the upper bounds 


ttsq < 3.70925 [98| and a™ < 3.55449 |99| have also been established, the value of a™ has 


now been rather well restricted. 


5.3 a and a for Strips of the Triangular Lattice 

For a strip of the triangular lattice, formed by starting with a square strip and adding 
diagonal edges, say from the upper left to the lower right vertices of each square, with free 
longitudinal boundary conditions, an elementary calculation yields 

a{tri, Ly = 2,Lx = m, FBCy, FBCx) = 2 ■ 3^"* . (5.3.1) 

For cyclic and Mobius strips of the triangular lattice with width Ly = 2, from the chromatic 


polynomials in |^, we hnd 
a{tri, Ly = 2,Lx = m, FBCy, PBCx) = —2 


7-h /7- a/Ts y 

2 ) 2 ) 


+ 5 (5.3.2) 


a{tri, Ly = 2, Lx = m, FBCy, T PBCx) = 3 "^”^ -|- 


vTS . 


7 a/Ts /7 - \/l3 

^r-) n^^) 


Using any of (|5.3.1|) , (|5.3.2|) , or (^.3.3|) , we compute 

a{tri, Ly = 2, FBCy) = 3 


(5.3.3) 

(5.3.4) 


For Ly = 3, using we get 

a{tri, Ly = 3, Lx = m, FBCy, PBCx) = —13 -|- 5 


3”^ + 


9 + 




(5.3.5) 


with 


Rm — ^ ^ ( ^t3c,j 
i=i 


(5.3.6) 
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(5.3.7) 


where the Xtscj, j = 1, 2, 3, are the roots of the equation 

+ 33^2 + 201^ + 324 = 0 

whence 

/43 + V1417 \h3 

a{tri, Ly = 3, FBCy) = - j = 3.4290909... (5.3.8) 

and, for Ly = 4, 5, 

a{tri, Ly = 4, FBCy) = {\tr^xyA,maxf^^ = 3.665350... (5.3.9) 

where Xtrixy 4 ,max IS the largest root of the quartic equation 

^4 _ 217^3 ^ 5950^2 _ 57953^ + 186624 = 0 (5.3.10) 

a{tri, Ly = 5, FBCy) = ( 

^trixyb^max) ^ 3.81486986... (5.3.11) 

where Xtxy 5 ,max is the root with the largest magnitude of the equation 

+ 1160^® + 330071^^ + 39854484^® + 2509883184^'^ + 87798785472^^ + 1684802267136^'^ 

+ 16500631191552^2 ^ 73085995450368^ + 104764180267008 = 0 


(5.3.12) 


For strips of the triangular lattice with cylindrical boundary conditions, one has 


a{tri, Ly = 3,Lx = m, FBCy, FBCx) = 6 ■ ( 71 )™" ( 5 . 3 . 13 ) 

a{tri, Ly = 3 , FBCy) = ( 71 )^/^ = 4 . 1408177 ... ( 5 . 3 . 14 ) 

For Ly = 4 from the chromatic polynomial calculated in we get 

a{tri, Ly = 4 , FBCy) = ( 352 )^/^ = 4 . 3314735 ... ( 5 . 3 . 15 ) 

Further, using the chromatic polynomials calculated in [|^, we have 

a{tri, Ly = 5 , FBCy) = (897 + 339^5 )^/'^ = 4 . 403125 ... ( 5 . 3 . 16 ) 

a{tri, Ly = 6 , FBCy) = {Xtxpy6,max)^^'' = 4 . 435287 ... ( 5 . 3 . 17 ) 

where Xtxpy6,max is the largest root of the equation 


- 8789 ^^ + 9102104 ^'* - 1119801408^2 ^ 42084913152 ^ - 364032294912 = 0 . ( 5 . 3 . 18 ) 


The monotonicity of a{tri, Ly, FBCy) and a{tri, Ly, FBCy) as a function of Ly is evident in 
the table. One can, of course, calculate these quantities for larger values of Ly-, however, it 
is interesting to observe that our exact calculation of a{tri, Ly = 6, FBCy) is already within 
about 1 % of the exact value for the Ly —>■ oo limit, i.e., for the inhnite triangular lattice. 


which we have obtained via an evaluation of the W function in 43 


a{tri) = 4.474647... (5.3.19) 

This shows that the values of a{tri, Ly, FBCy) converge rather rapidly to the inhnite-Lj^ 
value as Ly increases. 
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Table 3: Values of adG}) for the infinite-length limits of strip graphs of the triangular lattice with width 
Ly vertices and free (F) or periodic (P) transverse boundary conditions. 


Ly 

BGy 

a({G} 

2 

F 

3 

3 

F 

3.429 

4 

F 

3.665 

5 

F 

3.815 

3 

P 

4.141 

4 

P 

4.331 

5 

P 

4.403 

6 

P 

4.435 

CXO 

P 

4.475 


5.4 a and a for Cyclic Clan Graphs 

From eqs. (|1.8.1|) and ( p..8.2|) we have 

aiG[iKrUjn]) = {-lY^P{G[iKrU jn], q =-1) 


= i-iy 


K-l - r),,,]" + = -1))' 


d=l 


|(2r)„]’" + 25 ^(-l)''|r,„( 2 r),,_„]’ 


d=l 


where we have used 
and 


^P,r,d{q — 1 ) — ( 1 ) (r—d) 


(5.4.1) 

(5.4.2) 

(5.4.3) 


lid (4 ^-1) ^ (-1)'^2 ford>l 

(the d = 0 coefficient being a constant, /to = !)• 

Since the point g = — 1 is in region Ri, the term Xp^r,d with d = 0 is dominant, and hence 

a(Gl{(Kr)},jn] = |(2r)M|‘/' . (5.4.4) 

For large r —>■ oo, this has the leading asymptotic behavior a{G[{{Kr)}, j"^]) ~ 4e“^r. 

6 Number of Spanning Trees for Two G[{Kr)rm L] Fam¬ 
ilies 
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6.1 General 


In this section we shall again consider only connected graphs G = (V, E) and for notational 
simplicity, we let 

n=\V\. (6.1.1) 

We denote the number of spanning trees of such a graph as Nst{G). Here we shall prove 
theorems that determine the the numbers of spanning trees of two G[{Kr)m, L] families of 
graphs, namely the cyclic clan graph family G[{Kr)m, L = jn] and the family G[{Kr)mi L = 
id]. In general Nst{G) is given by ( [1.5.1|) ; however, we shall actually use an alternative 
method of calculation based on the Laplacian matrix (e.g., JT^). We recall that the adjacency 
matrix A(G) of a graph G is the n x n matrix whose jfc’th element is the number of edges 
connecting vertex j with vertex k in V. Next, we recall 

Def. The Laplacian matrix Q(G) of this graph is given by 


Q(G) = A(G) - A(G) 


( 6 . 1 . 2 ) 


where A is the nxn diagonal matrix whose j’th diagonal entry is equal to the degree of the 
j’th vertex, Ajj = A{vj) and whose other entries are zero. 

Since the sum of the elements in each row (or column) of Q(G') vanishes, one of the 
eigenvalues of Q(G) is zero. Denote the remaining n — 1 eigenvalues by Ai,..., A„_i. Then a 

T|, [Tool ) 


basic theorem is (e.g.. 


Nst{G) = any cofactor of Q(G) 


1 


n 


n—1 

i=i 


(6.1.3) 

(6.1.4) 


For the recursive families Gm of graphs considered here, Nst{G) grows exponentially for 
large m, and hence a quantity of interest is the growth rate. In previous works calculating 
z{{G}) and Nst{G) for the n —^ cx) limits of regular lattice graphs ||101|| - [|T0^ the convention 
was used of defining this growth rate in terms of In and we shall follow this convention 
here: 


Def. Let G be a connected graph. Then the quantity 2 ;({G}) is defined by 

exp(z({G}))= hm[A 5 T(G)]i/W 


(6.1.5) 


Since we shall compare the expressions that we derive for Nst on G[{Kr)m, and 
G[{Kr)m,id] with upper bounds, we recall the statements of these bounds. A general upper 
bound is ||104|| 

1 / 2\E\ \ 

A^5t(G)<-(^) . (6.1.6) 

n\n — 1/ 

For a A-regular graph G, using the relation A = 2\E\/n, this implies the upper bound 

1 / nA ^ 


Nst{G) < — ( 
n V 


n 


(6.1.7) 
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and hence 


exp( 2 ({G})) < A . (6.1.8) 

A stronger npper bonnd for A-regnlar graphs with vertex degree A > 3 is ||105| , |106|| : 

2 Inn 


where 


Nst{G) < 


C'a = 


nA In A 
( A -1)^-1 


(CaT: 


[A(A - 2)]^/2-i • 
yielding the npper bonnd for a A-regnlar graph with A > 3 


exp(^({G'})) < C, 


A • 


(6.1.9) 


( 6 . 1 . 10 ) 


( 6 . 1 . 11 ) 


By expanding Ca as A —> cx), one sees that in this limit the npper bonnd (|6.1.11|) approaches 

O- 

For the comparison of a given growth rate of a A-regnlar graph with the two npper 
bonnds (|6.1.8|) and (|6.1.11|) , labelled as u.b.,j with j = 1, 2, we dehne the ratios 


R,{{G}) = 


exp(^({G'})) 

exp(^({G})„.b,j) 


( 6 . 1 . 12 ) 


6.2 Number of Spanning Trees in G[{Kr)rm jn] 

Here we shall calcnlate the nnmber of spanning trees for the graph G[(iFr)m, jhi]- Let a(j, j') 
be the r x r adjacency matrix between the vertices of and Kr^,- The non-zero entries in 
this matrix are 

a(j,j)=Jr-Ir forl<j<m (6.2.1) 

and 

a(j,j + 1) = a(j + l,j) = forl<j<m (6.2.2) 

where J,. is the r x r matrix with all elements eqnal to nnity, and T is the r x r identity 
matrix. The Laplacian matrix of the graph G[{Kr)m, j"^] is then 


Q{G[{Kr)mJn\) = [(3r - 1)T - a(j, j)] Im - a(j, j -M) ® - a(j 1, j) ® 

= (3rT - J^) 0 0 R„ - (g) R^ (6.2.3) 

where Rm is the m x m matrix 


/ 0 1 0 
0 0 1 


0 0 0 
\ 1 0 0 


0 \ 
0 


1 

0 / 


(6.2.4) 
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XiQiG[iKrUjn])) = 


Lemma 3 The eigenvalues of Q{G[{Kr)m, jn]) are 

3r — (r + = 2r[l — cos(27rj/m)] for j = 0,1, • • 

3r with multiplicity (r — l)m. 

(6.2.5) 

Proof We hrst observe that Rm can be diagonalized by the similarity transformation 
SmRmS”^ generated by the matrix Sm with elements 

(S,„)jfc = for j,/c = 0,1,-1, (6.2.6) 

where * denotes the complex conjugate, and i = Therefore, Rm has the eigenvalues 


Aj(R,„) = for j = 0,1, ■ ■ •, m - 1. 

The result in ( |6.2.5|) then follows. □ 

Theorem 5 

NsT{G[iKrUjn]) = . 

Proof From (|6.1.4|) we obtain the result 


(6.2.7) 


( 6 . 2 . 8 ) 


1 

rm 

1 


m—1 


NsT{G[{Kr)mJn]) = —(3r)('’ JJ 2r[l - cos(27rj/m)] 

[r(3r)' 


i=i 

m—1 

\r—lim 


■ 


JJ 4sin2(j7r/ 
i=i 


(6.2.9) 


To evaluate the product, we use the relation 


m—1 


2 ^ TT sm(e + —) = sm(me) 

m 

j=0 


( 6 . 2 . 10 ) 


taking the limit e —>• 0 and applying L’Hospital’s rule. This yields NsT{G[{Kr)m, jn]) = 
(r^m)“^[r(3r)^“^]™'m^, which, in turn, gives the result in (|6.2.8|) . □ 

Corollary 5 Consider the number of spanning trees on the graph G[{Kr)m, j'^] ^-nd 
take the limit m ^ oo. The quantity measuring the growth rate in this limit is given by 


exp(2:(G[(iF,.)oo, j^^])) = 3 


i-J 




( 6 . 2 . 11 ) 


Proof This follows from (|6.2.8|) and the dehnition (|6.1.5|) . □ 

For the hrst nontrivial case, r = 2, eq. (|6.2.11|) yields exp{z{G[{K 2 )oo, = 2a/ 3 ~ 
3.46410, in agreement with our result in eq. (A. 16) of 


, m 
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Table 4: Values of eyi^ {z{G\ {Kr)oQ , jn])) , abbreviated as e^, from ( 6.2.11 ) for 2 < r < 10 and comparison 
with the upper bounds ( |6.1.8[) and ( |6.1.1l| ) via the ratios Rj, j = 1,2 given by ( 6.1.12 ). 


r 


Ri 

R2 

2 

3.464 

0.693 

0.786 

3 

6.240 

0.780 

0.838 

4 

9.118 

0.829 

0.871 

5 

12.041 

0.860 

0.894 

6 

14.988 

0.882 

0.910 

7 

17.950 

0.897 

0.921 

8 

20.920 

0.910 

0.931 

9 

23.897 

0.919 

0.938 

10 

26.879 

0.927 

0.944 


For large r, the ratios Ri, i = 1, 2 both approach unity: 


R,{G[{KrU,jn]) = 




1 + — Insjr ^ 

+ (^ - ^ln3 + ^(ln3)^)r-2 + 0{r-^) 

1 - 0.76528 r-^ + 0.34838 + 0{r-^) 


( 6 . 2 . 12 ) 


R2{G[{KrU,jn]) = 


1 + — ln3jr ^ 

+ (^ - ^ln3 + ^(ln3)2)r-2 + 0(r-3) 

1 - 0.59861 + 0.345835 + 0(r"3) 


(6.2.13) 


This tendency of the growth rates to approach the upper bounds as the vertex degree in¬ 
creases was also found in ||103|| for regular lattices. 

In Table ^ we list the growth rate ( |6.2.11|) for 2 < r < 10 and compare it with the upper 
bounds (|6.1.8|) and (|6.1.11|) . (We do not list the lowest case, r = 1, since Nst grows linearly 
rather than exponentially in m for this value of r.) 


6.3 Number of Spanning Trees in G[{Kj.)rm 'id] 

Here we calculate the number of spanning trees in the family of graphs G[{Kr)ni, id]. In |]^ 
and [^, 1^ studies were carried out of the chromatic polynomials of the family G[(iFj,)m, L] 
for general linkage L, and the case of the identity linkage L = id was called the bracelet 
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in 


graph, G[{Kr)m,,'id] = Parenthetically, we note that the chromatic polynomial for 

G[{Kr)m,id] was computed for r = 2 in 0, for r = 3 
sequently, by different methods, in 1^, and for r = 5, 6 
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m 


for r = 4 in |]^ and sub- 
|]. The graph G[{Kr)m,id] 
has \V\ = mr, \E\ = (l/2)mr(r + 1) and is a A-regular graph with uniform vertex degree 
A = r + 1. For our calculation, we hrst note that the dehnition of a(j, j) remains the same, 
as given in (|6.2.1|) and 


+ 1) = a(j + 1, j) = Ir for 1 < i < m 


(6.3.1) 


The Laplacian matrix for this the G[{Kr)mGd] graph is 

Q{G[{Kr)m,id]) = [(r + 1)1,. - a(j, j)] ® - a(j, j + 1) 0 R„ - a(j + l,j) 0 

[(r + 2)1,. - J^] ® - I,. (8) R,„ - I,. ® R^ (6.3.2) 


Lemma 4 The eigenvalues of Q{G[{Kr)m,id]) are 

r r + 2 - (r + ^ ^ 2[1 - cos(27rj7m)] 

X{Q{G[{KrUtd]) = \ 

[ r + 2 - ^ =r + 2[l- cos(27rj7m)] 

The proof proceeds in the same way as before, and is omitted. 

Theorem 6 The number of spanning trees in G[{Kr)mGd] is 


for j = 0,1, ■ • •, m — 1 

with multiplicity (r — 1). 
(6.3.3) 


NsT{G[iKr)m,td]) = 





w'^-2+uj-”^Y-1 


r—1 


for even m 
for odd m. 


(6.3.4) 


where 


u = 


r + 2 + v/r(r + 4) 


(6.3.5) 


Proof We hrst use the lemma (|6.3.3|) with (|6.1.4|) to obtain 

... rPR—1 1 m—1 


NsT{G[iKr)mGd]) = — 

rm 


— 


L j = l 

m—1 

r-2 ' 


Yl 2[1 - cos(27rj/m)] JJ 

L -I k=0 

m—1 

[r + 2 — 2 cos(27r/c/m)] 


r + 2[1 — cos(27rfc/m)] 


r—1 


'-k=l 


r—1 


(6.3.6) 


For even m, the product can be evaluated by using the identity 

— 1 ^2m _ 1 


7^ — 2a;cos(fc7r/m) + 1] = 


k=l 


— 1 


(6.3.7) 
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and setting a: = a; to get 


m—1 


pm/2—1 


[r + 2 — 2 cos(27r/c/m)] = (r + 4) [r + 2 — 2 cos(7r/c/(m/2))] 


k=l 


1 2 


= (r + 4) 


= (r + 4) 


k=l 


1_ m 

U! 2 


o; 


u 


2-1 


o; 


m/2 _ 


U! — U 


-1 


For odd m, we nse the eqnation 


— 2ci; cos{2nk/{2m + 1)) + 1] = 


u 


2m+l _ 


k=l 


Ci; — 1 


(6.3.8) 


(6.3.9) 


with the same snbstitntion (|6.3.5|) for cn, to obtain 


m—1 


[r + 2 — 2 cos(27r/c/m)] 


k=l 


(m-l)/2 

[r + 2 — 2 cos(27r/c/m)] 

k=l 


UJ 


UJ 


(l-m)/2 


W™- 1 


o; — 1 


2 -\- u] 


oj — 2 -\- u] ^ 

Combining eqs. ( |6.3.6|) , (|6.3.8|) and (|6.3.10D , we then have the result given in (|6.3.4|) . □ 


(6.3.10) 


We remark that the generating function of for even m is 1/(1 — (r + 2)(/ + f/^), 


and the generating function of -I /2 for odd m is (1 + C)/(l — (r" + 2)(/ + (/^). 


Corollary 6 Consider the number of spanning trees on the graph G[{Kr)miid] and 
take the limit m —>■ 00 . The quantity measuring the growth rate in this limit is given by 


,1- 


uj i.e., 


exp{z{G[{Kr)oc, id])) = 


r + 2 + 


n 1-2 


(6.3.11) 


Proof This follows from (|6.3.4|) and the dehnition (|6.1.5|) . □ 

For r = 1, exp( 2 :(G[(iFi)oo,/d])) = 1, reflecting the fact that Nst = m for G[{Ki)m,id] = 
Gm, which is subexponential growth in m. For r = 2, 

exp{z{G[{K2)oc, id])) = (2 + Vs )^^^ ~ 1.93185 (6.3.12) 

in agreement with our result in eq. (D.21) of |]^. For r = 3, 


/ 5 -I- \/2T \ 2/3 

exp(2;(G[(iC3)oo,/d])) = - - -j ~ 2.84207 (6.3.13) 
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Table 5: Values of exp(z(G 
with the upper bounds ( |6.1. 




(Kr)oo ,id\)), abbreviated as e^, from (|6.3.1l[) for 2 < r < 10 and comparison 
and ( |6.1.1l| ) via the ratios Rj, j = 1,2 given by ( 6.1.12 ). 


r 


Ri 

R 2 

2 

1.932 

0.644 

0.837 

3 

2.842 

0.711 

0.842 

4 

3.751 

0.750 

0.851 

5 

4.664 

0.777 

0.860 

6 

5.582 

0.797 

0.867 

7 

6.505 

0.813 

0.874 

8 

7.433 

0.826 

0.879 

9 

8.365 

0.836 

0.884 

10 

9.301 

0.846 

0.889 


in agreement with our result in eq. (A.88) of 


For large r, the ratios Ri, i = 1,2 both approach unity in the manner indicated below: 
Ri{G[{Kr)oo,id]) = l + (l-lnr)r“^ 


—4 — In r + - In r ) r + O ( 


/Inr \ 3 

r 


(6.3.14) 


R2{G[{Kr)oc,jn]) = l+(^^-lnrjr^ 

+ (-y - ^Inr+ + (6.3.15) 

In Table |] we list the growth rate (|6.3.11|) for 2 < r < 10 and compare it with the upper 
bounds (|6.1.8|) and ( |6.1.11| ). As expected, the approach of Ri and R 2 to unity as r increases 
is less rapid for G[{Kr)cx),id] than for G[{Kr)oo, jn] because of the fact that the family with 
L = id has a lower vertex degree, A = r + 1 than that of the family with L = jn, for which 
A = 3r — 1. 

Acknowledgments. R.S. thanks Profs. Dominic Welsh and Marc Noy for their organiza¬ 
tion of the stimulating CRM Workshop on Tutte Polynomials. This research was partially 
supported by the NSF grant PHY-0098527. 


7 Appendix 

In this appendix we list the cubic and quartic equations that yield the for T(G[(iV 3 )m, jn],x, y) 

that were not already given in the text. The AT, 3 ,i,j, j = 1, 2, 3,4 are solutions of the quartic 
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equation 


+ ^42'C^ + &41'C + &40 — 0 


(7.1) 


where 

643 = — (y® + 3y^ + 6y'^ + lOy® + 15y^ + 3a;|/^ + 21 ?/^ + Gxy"^ + 31y^ 

+3x‘^ + 9xy + 39y^ + 21x + 45y + 36) (7.2) 

642 = 3y{xy^^ + Axy^^ + y^^ + x'^y^ + lOxy^ + 5y^^ + 3x‘^y'^ + 23a;|/® 

+15|/® + 6x^y^ + 39xy'^ + 31y^ + lOx'^y^ + 53xy^ + A2y'^ + Ux^y* 

+&2xy^ + 41|/® + 3x^y‘^ + 17x^1/^ + Qbxy* + 28?/^ + 3x^y + 30x^1/^ 

+66x1/^ + 8|/"^ + 21x^1/ + 67x1/^ — 13y^ — 6x‘^ + Qxy — 38y‘^ — 30a; — 66y — 36) 

(7.3) 

&41 = —9y^{y + l){xy^^ + x^y^ + 5xy^^ + 2x^y'^ + Ax'^y^ + 16xy^ + 2x^y^ 

+7x^y’^ + 38xy^ + x^y^ + 7a;^?/® + 37xy'^ — 2y^ + 9a;^|/® + 72xy^ — lly’^ 

+lAx^y'^ + 81xy^ — 23y^ + 3x^y‘^ + 20 a;^|/^ + 733;!/"^ — Aly^ + 6 x^ 1 / + 2Qx^y‘^ 

+AAxy‘^ — Q3y'^ — 3x^ + 3x^y — 8Ay‘^ — 18x^ — A8xy — 85y'^ — 33x — 61y — 18) 

(7.4) 

640 = ‘^7y^{y + ifiy"^ + y + l){xy- l)(a;V + 2a;^l/^ + ‘2xy‘^ - x^y"^ 

—3xy^ — y‘^ — ^xy"^ — 3y^ + 2x + Ay + 2) . (7.5) 

The \T, 3 fi,j, j = 2,3,4, solutions of the cubic equation 

+ &32'C^ + &31'C + &30 = 0 (7.6) 

where 

632 = -{y^ + 3/ + 6|/^ + 10/ + I5y^ + 3xy^ + 21y^ + x^ + 6xy‘^ + 28y^ 

+9x‘^ + 16xy + 36y^ + 26a; + 38y + 24) (7.7) 
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hi = + 12xy^^ + Sx^y'^ + 9x‘^y^ + 36xy^ + 3y^^ 

+6x^y^ + 27x‘^y'^ + 90xy^ + 12y^ + lOx^y^ + hAx^y^ + lh3xy'^ + VJy^ 

+ I3x^y^ + 3Ax^y^ + 194:xy^ — Qy'^ + 3x'^y‘^ + 21x^y^ + 102x^1/^ + 13Axy^ 

—62y^ + 6x^y + AOx^y'^ + 9Qx^y^ + llTa;?/"^ — 13%® — 3a;^ + 18x^1/ + Ibx^y"^ 

+ lQxy^ — I92y^ — 2Ax^ — 3Qx^y — Sbxy"^ — 211y^ — 69x‘^ — 156xy — 18 %^ 

-84a; - 120y - 36) (7.8) 

ho = —^y'^iy + + 3x‘^y^ + 6x^y'^ + 3x'^y^ + 7x^y^ + 3x‘^y'’^ 

^x'^y'^ — 6x^y^ — Idx'^y^ — Axy'^ — 3x‘^y^ — Ibx^y'^ — 31x‘^y^ — 12xy^ 

—llx^y^ — lAx'^y'^ + 7xy^ + 5y^ + x'^y + 9x^y‘^ + 2lx‘^y^ + 33xy‘^ + 9y^ 

+6x^y + 23x‘^y‘^ + 23xy‘^ + y^ — 2a;^ — x^y — Axy'^ — lly^ — Qx'^ — 14:xy 
-12y^ -Ax- Ay) . (7.9) 
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